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ABSTRACT 



(N 

^ ! Motivated by the low-coUisionality of gas accreted onto black holes in Sgr A* 

O '■ a-nd other nearby galactic nuclei, we study a family of 2D advective accretion 



solutions with thermal conduction. While we only impose global inflow, the 



^, 

O , accretion flow spontaneously develops bipolar outflows. The role of conduction 

Q ! is key in providing the extra degree of freedom (latitudinal energy transport) 

necessary to launch these rotating thermal outflows. The sign of the Bernoulli 

Ph| constant does not discriminate between inflowing and outflowing regions. Our 

O I parameter survey covers mass outflow rates from ~ to 13% of the net inflow rate, 

'^ I outflow velocities from ~ to 11% of the local Keplerian velocity and outflow 

^ ' opening angles from ~ to 60 degs. As the magnitude of conduction is increased, 

outflows can adopt a conical geometry, pure inflow solutions emerge, and the 

rN I limit of 2D non-rotating Bondi-like solutions is eventually reached. These results 

c^ I confirm that radiatively-inefficient, hot accretion flows have a hydrodynamical 

propensity to generate bipolar thermal outflows. 

Subject headings: accretion, accretion disks - conduction - black hole physics - 
hydrodynamics 



1. Introduction 

Over the past decade. X-ray observations have made it increasingly clear that black 
holes are capable of accreting gas under a variety of flow configurations. In particular, 
there is convincing evidence that a distinct form of hot accretion exists at sub-Eddington 
rates, which contrasts with the classical "cold and thin" accretion disk scenario (Shakura 
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& Sunyaev 1973). Hot accretion appears to be common in the population of supermassive 
black holes in galactic nuclei and during the quiescent phases of accretion onto stellar-mass 
black holes in X-ray transients (e.g. Narayan et al. 1998a; Lasota et al. 1996; Di Matteo et 
al. 2000; Esin et al. 1997, 2001; Menou et al. 1999; see Narayan, Mahadevan & Quataert, 
1998b; Melia & Falcke 2001; Narayan 2003; Narayan & Quataert 2005 for reviews). 

Theories for the structure and properties of hot accretion flows remain controversial, 
however. Inspired by the work of Shapiro, Lightman & Eardley (1976), Narayan & Yi (1994; 
1995a,b) derived self-similar ADAF solutions which emphasized the important stabilizing 
role of radial heat advection (see also Ichimaru 1977; Rees et al. 1982; Abramowicz et al. 
1988). Subsequent analytical work on hot accretion flows has emphasized the possibility of 
outflows, motivated by a positive Bernoulli constant (Narayan & Yi 1994; ADIOS: Blandford 
& Begelman 1999), and the potential role of convection (CDAFs; Narayan et al. 2000, 
Quataert & Gruzinov 2000a). Hydro dynamical and MHD numerical simulations, on the 
other hand, have greatly contributed to the subject by highlighting important dynamical 
aspects of the problem which are not captured by idealized analytical models (e.g. Stone 
et al. 1999; Stone & Pringle 2001; Hawley et al. 2001; DeVilliers et al. 2003; Proga & 
Begelman 2003; McKinney & Gammie 2004). 

The purpose of this work is to consider the possibility that thermal conduction, which 
has been a largely neglected ingredient, could affect the global properties of hot accretion 
flows substantially. In §2, we use existing observational constraints on a few nearby galactic 
nuclei to argue that hot accretion is likely to proceed under weakly-coUisional conditions 
in these systems, thus implying that thermal conduction could be important. In §3, we 
extend the original ID self-similar solutions of Narayan & Yi (1994) to include a saturated 
form of thermal conduction and we study the effects of conduction on the flow structure and 
properties. In §4, we generalize the 2D solutions of Narayan & Yi (1995a) in the presence 
of thermal conduction and flnd that this extra degree of freedom allows the emergence of 
rotating thermal outflows. In §5, we discuss our results in the context of recently established 
observational trends for hot accretion flows and related outflows. In §6, we comment on sev- 
eral limitations of our work, potential additional consequences of strong thermal conduction 
in hot accretion flows as well as directions for future work. 



2. Observational Constraints and Collisionality 

Chandra observations provide tight constraints on the density and temperature of gas at 
or near the Bondi capture radius in Sgr A* and several other nearby galactic nuclei. These 
observational constraints have been used before to estimate the rate at which the gas is 
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captured by the black hole in these systems, following Bondi theory (e.g. Loewenstein et al. 
2001; Di Matteo et al. 2003; Narayan 2003). Here, we use these same constraints (taken 
from Loewenstein et al. 2001; Baganoff et al. 2003; Di Matteo et al. 2003; Ho, Terashima 
& Ulvestad 2003) to calculate mean free paths for the observed gas. It should be noted that 
these constraints are not all equally good, in the sense that Chandra has probably resolved 
the gas capture radius in Sgr A* and M87 but not in the other nuclei. 

Galactic nuclei and their associated gas properties on 1" scales are listed in Table 1: ni" 
is the gas number density, Tyi is the temperature, i?i is the 1" size-equivalent at the nucleus 
distance, /i is the mean free path for the observed 1" conditions and i?cap is the capture radius, 
inferred from the gas temperature and the black hole mass in each nucleus. From Spitzer's 
(1962) expressions for collision times, a standard expression for the thermal speed and a 
Coulomb logarithm InA ~ 20, one gets the simple mean free path scaling / ~ 10"^(T^/n) cm 
(Cowie & McKee 1997; valid for both electrons and protons in a 1-temperature gas). We 
have used this relation to calculate /i values in Table 1. 

The inferred mean free paths are in the few hundredths to few tenths of the observed 
1" scales {li/Ri values in Table 1). Assuming no change in the gas properties down to the 
inferred capture radius, values of /i/-Rcap ^0.1 are deduced in four of the six nuclei. One 
could also reformulate these scalings by stating that the 1" mean free paths are systematically 
> lO'^'^-Rs, where Rs is the Schwarzschild radius (the typical length-scale associated with 
the accreting black holes). These numbers suggest that accretion will proceed under weakly- 
collisional conditions in these systems. 

The weakly-collisional nature of ADAFs has been noted before (e.g. Mahadevan & 
Quataert 1997), in the sense of collision times longer than the gas inflow time, but these are 
model-dependent statements. Direct observational constraints on the gas properties near the 
capture radius favor a weakly-collisional regime of hot accretion, more or less independently 
of the exact hot flow structure. Whether the gas adopts an ADAF, CDAF or ADIOS 
type conflguration once it crosses the capture radius, it is expected to become even more 
weakly-collisional as it approaches the black hole, since the relative mean free path scales as 
l/R oc i?~^/^~P in these flow solutions with a virial temperature proflle and a density proflle 
p oc i?-3/2+p. 

Tangled magnetic flelds in the turbulent accretion flow would likely reduce the effective 
mean free paths of particles. The magnitude of this reduction, which will depend on the 
fleld geometry, is currently unknown. We return to this important issue in our discussion 
section. Even if the effective mean free paths of electrons are reduced to values lea/ R ~ 1, 
thermal conduction may be important . It is therefore of interest to investigate its effects on 
the structure and properties of hot accretion flows in general. We begin with a study of ID 



-4- 



ADAF models with conduction in §3 and then move on to study a 2D version of these hot 
accretion flows in 54. 



3. One-Dimensional Hot Accretion Flows with Saturated Conduction 

3.1. The Self-Similar Equations in One Dimension 

The standard Spitzer formula for thermal conduction applies only to gas well into the 
collisional ("fluid") regime, with a mean free path / ^ L, for any relevant flow scale L. 
We expect radial temperature variations on scales ~ R from self-similarity. In a weakly- 
collisional regime with / S> -R or even I r^ R, a. saturated (or equiv. "flux-limited") thermal 
conduction formalism should be adopted if one is to avoid unphysically large heat fluxes. We 
adopt the formulation of Cowie and McKee (1977) and write the saturated conduction flux as 
Fs = 5<l>sPCg, where <l>s is the saturation constant (presumably ~ 1), p is the gas mass density 
and Cs is its isothermal sound speed. With this prescription, the largest achievable conduction 
flux is approximated as the product of the thermal energy density in electrons times their 
characteristic thermal velocity (assuming a thermal distribution and equal electron and ion 
temperatures; see Cowie and McKee 1977 for details). Because it is a saturated flux, it 
no longer explicitly depends on the magnitude of the temperature gradient but only on the 
direction of this gradient. Heat will flow outward in a ID hot accretion flow with a near-virial 
temperature proflle, hence the positive sign adopted for Fg. 

With this formulation for conduction, the steady-state 1-temperature self-similar ADAF 
solution of Narayan & Yi (1994) can be generalized to include the divergence of the saturated 
conduction flux in the entropy-energy equation (see also Honma 1996 and Manmoto et al. 
2000 for investigations of the role of "turbulent" heat transport in one-dimensional ADAF- 
like flows). Adopting the same cylindrical geometry and notation as Narayan & Yi (1994), 
we consider a "visco-turbulent" , differentially rotating hot flow which satisfles the following 
height-integrated equations for the conservation of momentum and energy 

'"" dR ~ pRHdR V ^ dRJ ' ^ ^ 

^^P^'^lR - ^ ^K KdR) -^dR^^^'^' ^^^ 

supplemented by the continuity equation, which is rewritten as M = —47iRHvj.p. The 
last term in equation (3) is the additional divergence of the saturated conduction flux. In 
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these equations, R is the cylindrical radius, VLk is the Keplerian angular velocity around 
the central accreting mass, v^ is the gas radial inflow speed, Vt is its angular velocity, p 
is the mass density, Cg is the isothermal sound speed, s is the gas specific entropy, T is its 
temperature, H = Rcs/vk is the fiow vertical thickness and / < 1 is the advection parameter 
(/ = 1 for negligible gas cooling). A Shakura-Sunyaev prescription has been used to capture 
the "visco-turbulent" nature of the flow, with an equivalent kinematic viscosity coefficient 
u = ac^./Qx, where a is the traditional viscosity parameter. 

We look for solutions to these equations of the form 

p = po/2~'/', Vr = VqR-^'\ n = fio^"'/', cl = c%R-\ (4) 

and use the notation vk = VkoR~^^'^ = R^k = R^koR~^^^- We find that the three un- 
knowns Vq, Qq and c^o satisfy the relations 

^vl - Sl?,-„ + SIS + |c;„ = 0, (5) 

.0 - -^f . (6) 

9a2 4/5 5/3-7 \ 2 2 40$s 3 

-Vko = 0, (7) 

while the density scale, po, is fixed by the accretion rate, M. 

Contrary to the original ADAF solution, the energy relation (Eq. [7]) is not a quadratic 
in {cso/vkoY because of the extra saturated conduction term. We solve this fourth order 
polynomial with a standard numerical technique (Press et al. 1992) and discard the two 
imaginary roots as well as the real root with negative Cgo, which is unphysical because of a 
conduction fiux going up the temperature gradient. 



3.2. Results: Hot Accretion in One Dimension 

The main features of the new self-similar solutions with saturated conduction are il- 
lustrated in Figure 1, as a function of the value of the saturation constant, $s. We are 
showing results for two distinct solutions here, with 7 = 1.5 and 7 = 1.1, while we have fixed 
/ = 1 and a = 0.2 in both cases. Even though a value of 7 < 4/3 may not be physically 
meaningful, it is useful to illustrate more clearly how solutions depend on 7. 

Original ID ADAF solutions are recovered at small $s values. As $s approaches unity, 
however, the solutions deviate substantially from the standard ADAF. The sound speed, Cg, 
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and the radial inflow speed, —Vr-, both increase with the magnitude of conduction, while 
the squared angular velocity, Vt^, decreases. Formal solutions still exist when Vt^ becomes 
negative but they are obviously physically unacceptable. The solution with 7 = 1.5 reaches 
the non-rotating limit at log<l>s — ~2, while the solution with 7 = 1.1 reaches this limit for 
log$s — —0.8. Both solutions have sub-virial temperatures (c^ < Vk) and subsonic radial 
inflow speeds (— fr < c^), which appear to be rather general properties. Supersonic solutions 
may exist for large a values when 7 -^ 5/3. 

As the level of saturated conduction is increased, more and more heat flows from the 
hotter, inner regions, resulting in a local increase of the gas temperature relative to the orig- 
inal ADAF solution. Simultaneously, the gas adjusts its angular velocity (which reduces the 
level of viscous dissipation) and increases its inflow speed to conserve its momentum balance 
(increasing at the same time the level of advection). One can show from the energy equation 
that Vl? (X {1 — lo'cond/q'advl), whcre gcond and gadv are the magnitudes of the conduction and 
advection terms: solutions cease to exist once advection balances exactly the heating due to 
saturated conduction. We have found that variations in the advection parameter, /, have 
relatively minor effects on the solutions, as long as / is not ^ 1. The breakdown of the 
solutions (when Vl? -^ 0) occurs at lower values of the saturation constant <l>s for smaller 
a and larger 7 values. Solutions with Vl? = Q can be viewed as Bondi-like solutions with 
saturated conduction. 

Clearly, conduction can signiflcantly affect the structure and properties of hot accretion 
flows, even at a level well below the saturated flux value (i.e. for <l>s <^ 1). By analogy with 
the reasoning of Narayan & Yi (1995a) on convection in 2D self-similar ADAF solutions, one 
might expect conduction to be more important in the polar regions of hot accretion flows 
than at their equator. To better understand the role of conduction, it is therefore important 
to study it in a more realistic two-dimensional setting, without height integration. 



4. Two-Dimensional Hot Accretion Flows with Thermal Conduction 

4.1. The Self-Similar Equations in Two Dimensions 

We now extend the analysis of the previous section to two-dimensional hot accretion 
flows, following closely the methods of Narayan & Yi (1995a). We also discuss in Appendix 
C the relation between our work and that of Anderson (1987). Adopting spherical polar 
coordinates {r9(f)), we consider axisymmetric, steady-state flows {d/dcp = d/dt = 0). We 
deflne Keplerian angular and linear velocities as, respectively. 



Qk = {GMY'^r~'"\ 



'■(Pv) = 


= 0, 




^ ■ V)v = 


- —VP-nir + ^V-a, 
P ik 


(TVs) - 




l(^.^)_|(V.v)2 
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t;K = f^K^, (9) 

and we use the same a-prescription as before to describe the visco-turbulent nature of the 
flow: 

z/ = acl/VL^, (10) 

where c^ is the isothermal sound speed (c^ = Pj p\ 

The fluid equations obeyed by the visco-turbulent accretion flow are written as follows: 

(11) 
(12) 

+ -V-(AVT), (13) 

P 

where a is the viscous stress tensor (divided by the dynamic viscosity, pv). The notation 
for the other flow variables is the same as in §3. It should become clear, later on, that our 
main results are not expected to depend critically on the exact form of the viscous stress 
tensor adopted. To allow a direct comparison with our results in §3 and for consistency 
with Narayan & Yi (1995a), we use the same advection parameter / to quantify the cooling 
efficiency of the flow. We will see, however, that this parameter is in fact ill-defined for 
solutions with conduction. Fully developed forms of these equations (with a few typos and 
no thermal conduction term) can be found in Narayan & Yi (1995a). 

Contrary to the formulation adopted in §3, we do not write the conduction fiux in a 
saturated form here, but return instead to a traditional form where the fiux depends linearly 
on the local temperature gradient. Justifications for this choice are provided in Appendix A. 
We adopt a form of the thermal conductivity coefficient which preserves the radial self- 
similarity of the solutions: 

A(r) = Aor-i. (14) 

Following Narayan & Yi (1995a), we consider separable solutions satisfying 

Vr{r,e) = VK{r)vr{e), (15) 

veir,e) = 0, (16) 

v<p{r,e) = VK{r)n{e)sme, (17) 

p{r,e) = p{e)T-'/\ (18) 

Cs{r.e) = VK{r)cs{e). (19) 

In Narayan & Yi (1995a), the sin^ in the equivalent of Eq. (17) is missing, in what appears 
to be a typo. Note that we make the same f e = assumption as Narayan & Yi (1995a). 



With these forms, solutions automatically satisfy the continuity equation, Eq. (11). The 
problem is reduced to a one-dimensional differential system along the polar coordinate, 9. 
We are left to solve the following four coupled differential equations for the dimensionless 
variables Vr{0),^{0), p{0), and Cs{6): 



,2 /in2 



P[ --<-sin"^f]^ 



sin 9 cos 9pQ'^ 



- sin 9pVr^ 






dP / dP 3 dvr , 



a 



3„^ ^r.dn ^ d f dVL 

--Pn + 3cos9P— + sm9— ( P— 

4 d9 d9 \ d9 



(20) 
(21) 
(22) 



3 , 

--e Vr 
2 



a{3vl + ( —^ I + sin^ 



+K^ 



?n=+(^V 

4 \d9 



2 1 d f . ^dc^ 
'^ + ^9d9[''''^l9 



(23) 



For convenience, the parameters 






15/3-7 



/ / 7-1 



K = T-^O, 



(24) 
(25) 



were introduced in the energy equation (Eq. [23]), where 7 = Cp/c^ is the adiabatic index of 
the gas. 

This is a seventh-order differential system of four variables. The three momentum 
equations (Eqs [20-22]) are identical to those presented in Narayan & Yi (1995a). The 
energy equation (Eq. [23]) differs in that extra terms corresponding to the radial and polar 
contributions to heat conduction are now present. This qualitative difference turns out to 
be crucial in providing the extra degree of energetic freedom necessary for inflow or outflow 
at the pole. While the original energy equation seemed to require that rotating solutions 
be static {vr = 0) along the polar axis because of the positivity of all the viscous terms (no 
rotating solutions with f ^ 7^ at ^ = were found by Narayan & Yi 1995a), the addition 
of thermal conduction terms allows for radial velocities at the pole which are generally non- 
zero. Furthermore, for outflows to be energetically favored, the polar contribution to heat 
conduction (second term in the last square bracket in Eq. [23]) must be negative enough to 
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overcome the positive sum of the radial conduction term and all the viscous terms on the 
RHS of Eq. (23). When that happens, advection can act as an effective heat source, with a 
positive sign for the radial velocity, Vr (on the LHS of Eq. [23]). As we shall see below, this 
sign change for Vr preferentially happens in polar regions of the flow and only for certain 
combinations of model parameters. 

Boundary conditions for the system of differential equations are obtained by requiring 
that all flow variables be symmetric and continuously differentiable across the pole and the 
equator. This corresponds to the following conditions at ^ = and 6 = 7r/2: 

de de de de ' ^ ' 

As noticed by Narayan & Yi (1995a), these conditions are not all independent from each 
other, given Eqs. (20-23), so that the system ends up not being over-constrained. 

The net mass accretion rate, which one obtains from the volume integral of the continuity 
equation (Eq. [11]), is 

m = - I 2np{9)v{9) sin 9d9, (27) 



where the radial dependence has been omitted because of self-similarity. This integral sup- 
plies two additional boundary conditions to our system because we enforce it by solving the 
corresponding differential equation for dm/d9 and imposing the values of m at ^ = and 
9 = 7r/2. The m integral then uniquely determines the scaling of p{9) in our solutions. In 
what follows, we have used m = 1 (from ^ = to 7r/2) consistently in all the solutions. 

We find numerical solutions with an implementation of the relaxation code described in 
Hameury et al. (1998). We do not implement the adaptive mesh feature but simply solve 
the set of equations on a grid of 400-1000 points equally spaced in 9. After recovering the 
original 2D ADAF solutions of Narayan & Yi (1995a), we explored the parameter space of 
solutions with conduction in detail, in terms of the model parameters a, e' and Aq. 



4.2. Results: Hot Accretion in Two Dimensions 

4-2.1. Dynamical Structure of the Flow 

Figure 2 shows latitudinal profiles (^ = at the pole) of the four dynamical variables 
in our model, for solutions with a = 0.1, e' = 10 and several values of Aq. Except for the 
density, p (which is normalized to m = 1), the radial velocity, Vr, the angular velocity, Q, 
and the isothermal sound speed squared, cj., are all scaled to the local Keplerian values. For 
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small enough values of Aq (=10~^), solutions remain very close to the original 2D ADAF 
solutions of Narayan & Yi (1995a), with Vr (nearly) zero at the pole. Note that this solution, 
with e' = 10, is denser and cooler at the equator than at the pole and it corresponds to 
something like a "thin disk" surrounded by a hot coronal atmosphere. 

As the value of Aq is increased, the solutions start deviating substantially from the orig- 
inal ADAF, with faster radial inflow at the equator and signiflcant radial outflow appearing 
in polar regions. The solution with Aq = 1 in Fig. 2 illustrates this transition clearly. The 
flow also becomes hotter (except around the poles), is rotating more slowly (being more 
pressure supported) and has a more uniform density proflle, with a reduced density at the 
equator: it is gradually approaching spherical symmetry. As we shall see below, some of 
these variations in properties are shared by the original ADAF solutions if the flow is forced 
to be more advective (by reducing the value of e'). 

As the value of Aq is increased further, another qualitative change in the solutions 
occurs, with inflow at the pole, outflowing regions away from the pole and still faster inflow 
at the equator. In other words, the outflow adopts a conical geometry. The solution with 
Aq = 1.6 in Fig. 2 illustrates this regime. For even larger values of Aq, the solution switches 
to a global inflow pattern (Aq = 2.0), with increasingly faster inflow in the polar regions, 
which ultimately overcomes the inflow speed at the equator (Aq = 2.35). At the end of this 
sequence (Aq ^ 2.35), the flow approaches the non-rotating, uniform density and uniform 
temperature limit: the 2D solutions with conduction become Bondi-like, in a close analogy 
with what we have seen for ID solutions with conduction in §3. In this limit, the value of 
p ~ 5rn/{6na) = 2.65 (see Appendix B). We note that along this entire sequence of solutions, 
maximum inflow and outflow speed remained everywhere subsonic. As we shall see below 
in more detail, most of these solution properties can be understood in terms of the relative 
contributions of radial and polar conduction terms to the flow energetic balance. 

Figure 3 shows a comparable sequence of latitudinal proflles in solutions with a reduced 
value of e' (= 2) and the same a = 0.1 as before. For small enough values of Aq (= 10~^), 
the solutions are again close to the original 2D ADAF solutions of Narayan & Yi (1995a), 
with a more spherically symmetric proflle than in the e' = 10 case, but still a Vr value which 
is (nearly) zero at the pole. As the value of Aq is increased, a sequence similar to that 
in Fig. 2 obtains, with a few important differences: inflow and outflow velocities, outflow 
opening angles and the critical Aq value at which the non-rotating limit is reached are all 
smaller than in the e' = 10 case. In addition, there is no conical outflow regime along the Aq 
sequence in this subclass of solutions. 

Figure 5 illustrates in a different way the flow geometry and the relative importance 
of the rates of mass inflow and outflow, for four of the solutions shown in Fig. 2 {a = 0.1, 
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e' = 10). Isodensity contours are shown by solid lines, while color contours trace the flux 
of mass per unit solid angle (= pv^). The scale of the cylindrical coordinates {Rz) used to 
label the axes, as well as the isodensity contour levels, are arbitrary because of the radially 
self-similar nature of the solutions considered. 

Figure 6 summarizes the most significant results for the existence and the properties of 
outflows that emerged from our systematic survey of the parameter space of hot accretion 
flow solutions with conduction. Focusing on the specific a = 0.1 case, the various curves 
in Fig. 6 cover the (e', Aq) space. The four panels show (a) the flow radial velocity at the 
pole, Vr{6 = 0), in units of the Keplerian value, (b) the angular velocity at the equator, 
Q{6 = 7r/2), in units of the Keplerian value, (c) the opening angle of outflowing regions, 
including conical outflow cases, and (d) the ratio of outflowing to net inflowing mass rate 
(which is also the accretion rate defined in Eq. [27]). 

A number of clear trends emerge from this parameter space survey. Solutions with small 
e' values (^1) do not exhibit outflows, but only pure inflow regimes. Conversely, solutions 
with the largest values of e' are the most efficient at producing outflows. Furthermore, even 
among the solutions with e' > 1, there is only a finite range of Aq values for which an outfiow 
regime exists. The larger e' is, the larger are the velocity, opening angle and mass fiux of the 
outfiow, and so is the maximum value of Ag/e' beyond which the outfiow ceases to exist. We 
also find that solutions with conical outfiows exist only for sufficiently high e' values. Finally, 
for large e' solutions, the outfiowing regime persists even to very small values of Ap/e' but 
outfiows become increasingly weaker, slower and narrower. It should not be a surprise that 
these various fiow properties scale with the ratio Ag/e', given the form of the energy equation 
(Eq. [23]). All these trends are consistent with the detailed results shown in Figs. 2 and 3. 

The increased efficiency of large e' solutions at producing outfiows can be interpreted 
as resulting from the larger latitudinal temperature gradients present in these solutions. 
This makes the contribution of the polar conduction term more important, relative to the 
radial conduction term (which is tied to a fixed radial temperature gradient set by self- 
similarity). Mass outfiow rates ~ 6% of the net infiow rate, outfiow velocities ^ 2% of the 
local Keplerian velocity and outfiow opening angles ~ 50 degs are reached in solutions with 
e' = 10 and Aq ~ 0.1. For even larger values of e', outfiow properties continue to follow the 
trends shown in Fig. 6. For instance, we have found e' = 50 solutions with mass outfiow 
rates < 13% of the net infiow rate, outfiow velocities < 3% of the local Keplerian velocity 
and outfiow opening angles ^ 60 degs. Numerical convergence becomes increasingly difficult 
at large e' values, however, and we have not explored this regime as systematically as the 
low e' one. 

Two additional features in Fig. 6 are worth emphasizing. All solutions reach the non- 
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rotating limit for values of Aq/c' ^ 0.24 and Vr ^ —0.06 f k. A derivation of these asymptotic 
values for the non- rotating limit is provided in Appendix B. We also find that there is a 
wedge in the Aq/c' vs. e' space, such that solutions do not exist for a growing range of Ap/e' 
as the value of e' is increased. Despite careful efforts to reach numerical convergence in these 
parts of the parameter space, we have failed to find solutions of sufficient numerical accuracy 
(or at all) in these regions. We have not clearly identified the origin of this difficulty. Let 
us stress that all the solutions shown in figures, and otherwise discussed in this work, have 
reached satisfactory numerical convergence. 

We find that the value of the viscosity parameter, a, affects quantitatively (but not 
qualitatively) the properties of hot accretion flow solutions with conduction. For larger val- 
ues of a, the radial inflow and outflow speeds are increased and the overall density scale is 
reduced, as expected from the original ADAF a scalings. Figure 4 shows a speciflc compar- 
ison of solutions with e' = 10, Aq/c' = 0.1 and a = 0.02, 0.1 and 0.2. The outflow speed 
is much larger in the a = 0.2 solution, reaching a value ~ 0.11 fx at the pole. Yet, the 
mass outflow rate (~ 9%m) is not much greater in the a = 0.2 solution than in the other 
two solutions (~ 6-7% m) because of the reduced density scale associated with faster radial 
velocities. The outflow opening angles are comparable in all three solutions. For consistency 
and conciseness, we focus on solutions with a = 0.1 below. 



4-2.2. Detailed Energy and Momentum Balance 

To understand better the various properties of these hot accretion flows with conduction, 
let us now focus on the detailed energy and momentum balance that they satisfy. 

Figure 7 shows latitudinal proflles of the advection term (a), on the LHS of the energy 
equation (Eq. [23]), and of the conduction (b) and viscous (c) terms, on the RHS of Eq. (23), 
for the same speciflc sequence of solutions as shown in Fig. 2 (a = 0.1, e' = 10). In the 
solution with negligibly small conduction (Aq = 10"^), the advection and the viscous terms 
balance each other exactly and heat advection acts to cool the flow (which corresponds 
to positive values in panel a). In the solution with Ag = 1.0, however, the magnitude of 
the conduction term is already very important, being comparable to or larger than the 
magnitude of the viscous term in most of the flow. In addition, conduction both cools the 
flow in the polar regions (which corresponds to negative values in panel b) and heats it 
up in the equatorial regions. As a result of this cooling of the polar regions (where the 
contribution from the viscous term is comparatively small, simply by geometry), the balance 
between these various terms permits a negative advection term in the polar regions, which 
corresponds to advective heating and outflow. 
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As the value of Aq reaches 1.6, the latitudinal profile of temperature starts to flatten 
significantly (see Fig. 2), thus reducing the contribution of the polar component of conduc- 
tion. The radial component, on the other hand, is strictly positive and can only increase in 
magnitude with the value of Aq (because of the radial self-similarity of the solutions). As 
a result, the net conduction term switches sign again in the polar regions (which are now 
heated by conduction), but it is able to remain somewhat negative at intermediate values of 
the polar angle, 9. This situation leads to the conical outflow regime previously identified. 

Finally, for even larger values of Aq(> 2), the radial component of conduction dominates 
everywhere over the polar component (which is further reduced by the continued flattening of 
the temperature profile; Fig. 2), thus forcing the net conduction term to act strictly as a heat 
source (as in the ID solutions of §3). The energy equation is then satisfied only if advection 
acts to cool the flow everywhere, which must then be globally inflowing. For the largest 
values of Aq considered, advection very nearly balances conduction as the viscous term drops 
precipitously and the non-rotating limit is reached. In good analogy with what we have 
seen for ID solutions with conduction in §3, the non-rotating Bondi-like limit corresponds 
to strict equality between conduction and advection. 

From this examination of energy balance in our radially self-similar solutions, it is clear 
that the key to thermal outflows is the existence of strong enough cooling in the flow for 
advection to be allowed to act as a heat source. This is preferentially achieved in polar 
regions, simply because these are regions where the rate of viscous dissipation must drop by 
geometrical constraint. 

Let us comment here on the fact that the advection parameter, /, is actually ill-defined 
in our solutions with conduction. Since / is associated with the viscous dissipation term 
in Eq. (23), it represents a form of cooling which is proportional to the rate of viscous 
dissipation. As a result, in solutions with an energy balance dominated by conduction (e.g. 
solutions with large Aq in Fig. 7), a small value of / (<^ 1) does not correspond to a globally 
efficient cooling anymore, as it used to in the original ADAF solutions. This undesirable 
situation would disappear in 2-temperature solutions with more realistic descriptions of 
cooling, such as the ones described by Narayan & Yi (1995b). Note also that this is not 
a critical flaw for our analysis: the energetics of all the solutions discussed here depends 
only on the value of e' = e//, so that large e' values can equally well be interpreted as 
corresponding to low 7 values. 

For completeness, we now examine in detail the momentum balance in these same 
solutions {a = 0.1, e' = 10). Figure 8 shows latitudinal profiles of various terms in the 
radial and polar components of the momentum equation (Eqs. [20] and [21], respectively). 
The magnitude of the various terms on the LHS and the RHS of these equations are labeled 
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with different colors in the panels of Fig. 8, for the same sequence of increasing Ag values 
as before. Note that we ignore the azimuthal momentum equation, which simply expresses 
angular momentum conservation. 

For a negligibly small level of conduction (Aq = 10^®), the r-momentum balance is 
dominated by pressure support against gravity in the polar regions, while rotational support 
becomes gradually dominant in the equatorial regions. As the value of Aq is increased, the 
flow becomes hotter, more spherically symmetric and the r-momentum balance tends to be 
more dominated by pressure support at all values of the polar angle, 9. The inertial and 
viscous terms are never very important for r-momentum balance along this Aq sequence. 
As for the ^-momentum equation, for all values of Aq considered, it is satisfied through a 
balance between latitudinal pressure gradients and projected centrifugal force, with little 
contribution from the viscous terms. Perhaps one of the most interesting property of the 
profiles shown in Fig. 8 is that they show explicitly that there is no distinction, from the point 
of view of momentum balance, between inflowing and outflowing regions of the flow. This 
confirms that outflows, when they exist, are purely determined by energetic considerations. 
This further justifies us classifying them as "thermal" in origin. 



4-2.3. The Bernoulli Parameter 

Much of the discussion in the literature on the likelihood of outflows from hot accretion 
flows is articulated around the value of the Bernoulli parameter, B^. (or, more generally, the 
Bernoulli function). Narayan & Yi (1994; 1995a) first emphasized the positivity of -Be in their 
solutions and raised the possibility that this would lead to bipolar outflows. Furthermore, 
Blandford & Begelman (1999) 's motivation and method to develop ADIOS is rooted in the 
notion that hot accretion flows with positive B^ must have outflows and that their structure 
and properties should be accordingly modified, leading to a mass accretion rate varying 
with radius. The relevance of the value of Be for the generation of outflows has also been 
discussed in the context of numerical simulations (e.g. Stone et al. 1999) and it has been 
put into question (e.g. Abramowicz et al. 2000). Since our hot accretion flow solutions 
with conduction produce outflows spontaneously, it is of considerable interest to study the 
relation between these outflows and the value of B^ in the solutions. 

Figure 9 shows latitudinal profiles of the dimensionless Bernoulli parameter, b{6) = 
Bc/{rVt]^'^ (same notation as Narayan & Yi 1995a) in a set of solutions with a = 0.1, 
7 = 1.5 and different values of / and Aq/c'. In addition to b{9), scaled latitudinal profiles 
of the radial velocity, Vr, are shown, for an easy identification of infiowing vs. outfiowing 
regions. The nine panels in Fig. 9 are organized in such a way that the three rows correspond 
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to e' ~ 10, 2 and 0.5, from top to bottom, with increasing values of Ag/e' from left to right. 

From this ensemble of solutions, one concludes that the positivity of h{9) does not in 
general discriminate between inflowing and outflowing regions in our solutions. Although 
the solution corresponding to / = 0.033 and Ag/e' = 0.1 would seem to indicate it does, this 
is in fact a coincidence since we flnd many counterexamples. Notice in particular in Fig. 9 
how solutions with h{9) changing from positive at the pole to negative at the equator can 
be alternatively (nearly) static, outflowing or inflowing at the pole. Similarly, solutions with 
globally positive values of h{9) can be globally inflowing without any difficulty. 

The additional panels in Figure 10 show two particularly striking examples of the loose 
relation between the sign of h{d) and the inflow/outflow properties of our solutions. In one 
solution, with Ag/e' = 0.16, the outflow has a conical geometry. While the value of h{d) 
decreases monotically with 9, the radial velocity is clearly not monotonic, with two sign 
changes. In the other solution, with Ag/e' = 0.02, the outflow is polar but there is a small 
regions at mid-latitudes where outflowing gas has negative h{9) value. Although this type 
of "violation" is generally marginal, it clearly demonstrates that there is no strict relation 
between the sign of h{9) and the inflow/outflow properties of gas in our solutions. It should 
be noted that for large enough values of Ag, one would expect the Bernoulli parameter (which 
is an adiabatic, inviscid quantity) to lose its physical significance, since the flow must then 
become strongly diabatic. We suspect that the h{9) < outflow region shown in the left 
panel of Fig. 10 and the globally inflowing solutions with globally positive b{9) values shown 
in Fig. 9 are all manifestations of this loss of signiflcance of b{9) in flows with signiflcant 
conduction. 



5. Relevance to Observational Trends 

It is difficult to evaluate, from self-similar solutions alone, what the radiative efficiency 
of hot accretion flows with conduction may be. In our solutions, the flow structure adjusts 
so that the divergence of the conduction flux generated in the innermost regions is balanced 
by an increased level of entropy advection. There will be no extra source of energy available 
in the presence of conduction (in fact less is available in our solutions, since the viscous 
dissipation rate drops with Q), but the dissipated energy is redistributed differently in flows 
with conduction than in flows without. This redistribution affects the local flow properties 
(density, temperature), which in turn affect its cooling properties. 

Using our ID solution scalings (§3) for simplicity, it is possible to estimate crudely the 
modifled level of emission expected relative to a standard ADAF. Free-free emission scales as 
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p^T^I'^ oc (csofg)""^ for a given value of the accretion rate, M. We find that this emission is 
effectively reduced, e.g. by as much as ~ 25% in the ID solution of Fig. 1 with 7 = 1.5 (and 
much more in the 7 = 1.1 solution). Even though conduction heats up the gas locally, the 
reduced density resulting from the larger inflow speed dominates, leading to a net decrease 
in the expected level of free-free emission. On the other hand, the very steep dependence 
of synchrotron emission on the electron temperature (e.g. Mahadevan 1997) suggests that 
hotter solutions (with conduction) may be radiatively more efficient from the point of view 
of that process. These are only rough scalings, however. In light of several complications 
expected in global solutions with conduction (see discussion in §6), it is unlikely that we 
can determine, from self-similar solutions alone, how the global radiative efficiency of a hot 
accretion flow may be affected by conduction. We postpone a thorough discussion of this 
important issue to future work and focus here instead on the possible relevance of the outflows 
found in our solutions to a variety of observationally established trends for accreting black 
holes. 

In recent years, much progress has been made in understanding the outflow /jet prop- 
erties of accretion flows around black holes in X-ray binaries and Active Galactic Nuclei. 
Various observational studies suggest that, in both classes of systems, there is a regime of 
steady outflow when black holes accrete well below the Eddington rate (for luminosities 
^ 0.01 — O.lLEdd), while the outflow regime becomes flaring/unsteady at higher accretion 
rates (e.g. Ho 2002; Gallo, Fender, Pooley 2003; Merloni, Heinz & Di Matteo 2003; Fender, 
Belloni & Gallo 2004; Nipoti, Blundel & Binney 2005; Sikora, Stawarz & Lasota 2006). Some 
of the latest results support the notion that black holes accreting at sub-Eddington rates 
do so through an advective form of accretion (e.g. Greene, Ho & Ulvestad 2005; Koerding, 
Fender & Migliari 2006), but there is some evidence in the "low/hard" state of at least two 
Galactic black hole systems (Miller et al. 2006) that a standard thin disk remains embedded 
in the hot flow (a situation which may favor outflows according to our large e' solutions). It 
is tempting to relate these various observational trends with our results on thermal outflows 
from steady, hot accretion flows with conduction. 

The association would be strengthened if one could show that thermal conduction is 
expected to be equally important in hot accretion flows around stellar-mass and supermassive 
black holes. After all, we made the case for weakly-collisional accretion conditions only for 
systems harboring supermassive black holes in §2. While the temperature proflle does not 
depend on black hole mass in self-similar solutions scaled in units of the Schwarzschild radius, 
density proflles do. To address this question, we compare the scalings for the coefficient of 
turbulent viscosity [y = aCg/fix) and the Spitzer coefficient of thermal diffusivity (x oc 
T^'^/p), using for simplicity the ADAF scalings in Schwarzschild units of Narayan et al. 
(1998b). The comparison of these two diffusivity coefficients is the relevant one since it 
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captures the relative magnitude of "viscous" and conduction effects on any scale of interest. 
We find that the two coefficients scale identically (linearly) with the mass of the central black 
hole, which suggests that conduction would indeed be equally important in hot accretion 
flows around stellar-mass and supermassive black holes. This conclusion is also consistent 
with the similar scalings with black hole mass of the accretion time and collision times 
discussed by Mahadevan &; Quataert (1997). It would be interesting to strengthen this 
comparison by using, as much as possible, observational constraints on hot accretion flows 
in X-ray binaries, as we have done for AGN in §2. 

Establishing on a firmer ground the relevance of the thermal outflows found in our 
solutions to observational trends for accreting black holes will clearly require additional 
work. In the meantime, we note that our solutions allow a direct scaling of the expected 
outflow power, per unit accreted mass, which is a quantity of obvious interest. Recalling that 
the Bernoulli constant, b{6), represents the asymptotic value of the energy of gas escaping 
from the central mass gravity, we can use it as a measure of the outflowing energy rate and 
compare that to the rate of energy release due to accretion. The outflow power is simply 
Pout = Montb{9)v^. The accretion power is Pace — O.lMinC^ , if we pick a radius of maximum 
energy release of ~ 3 Schwarzschild radii to flx the accretion efficiency to ~ 0.1 as is usually 
done. At 3Ps, v^ ~ c^/6. For a ratio of Mont/Mi^ ~ 0.05 (solutions with e' = 10 in Fig. 6) 
and a corresponding value of b{6) ~ 0.15 for outflowing regions (Fig. 9), we deduce a ratio 
Pout /Pace that reaches up to ~ 1% in our solutions. Therefore, provided that solutions with 
realistic cooling properties can be constructed with radiative efficiencies < 10~^ (i.e. less 
than a hundredth of the fiducial 0.1 value), it may be possible to explain the "jet-power 
dominance" of radiatively-inefficient accretion flows which has been recently emphasized by 
several authors (e.g. Nagar, Falcke & Wilson 2005; Koerding et al. 2006; Allen et al. 2006). 
Interestingly, values of Pout/-Paec ^ 1% as we just derived are comparable to the values of 
-PjetZ-pBondi of a few perceuts inferred by Allen et al. (2006). 

In this context, understanding how thermal outflows may relate to the relativistic jets 
of accreting black holes will likely be important. After all, outflow velocities decrease with 
radius as r~^/^ in our solutions (by virtue of self-similarity), which is a rather unappealing 
feature. Thermal outflows could be collimated and accelerated by their environment or, 
alternatively, they could contribute to the collimation and acceleration of narrower and 
faster jets, such as the ones emerging in relativistic MHD numerical simulations of adiabatic 
tori (e.g. DeVilliers et al. 2005; Hawley & Krolik 2006). 

One of the best systems available to us for the study of radiatively-inefficient accretion 
flows (and related outflows) is Sgr A*. In recent years, the focus of theoretical work to 
interpret the rich set of data has in fact shifted to include a jet or outflow component (e.g. 
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Falcke & Biermann 1999; Falcke & Markoff 2000; Melia & Falcke 2001; Markoff et al. 2001; 
Yuan, Markoff & Falcke 2002). Observations of the size (e.g. Bower et al. 2004; Slien et al. 
2005; Yuan, Slien, Huang 2006), variability (e.g. Yuan, Quataert & Narayan 2003; Yusef- 
Zadeli et al. 2006) and polarization (e.g. Agol 2000; Quataert & Gruzinov 2000b; Bower 
et al. 2003; Marrone et al. 2006) of this source all have tremendous potential to further 
constrain the properties of the hot accretion flow and its outflow, and they may be able to 
help us discriminate between the thermal outflows seen in our solutions and other outflow 
alternatives. Direct constraints on outflow geometries for accreting black holes are generally 
scarce, but it is interesting to note that Junor, Biretta & Livio (1999) observe an opening 
angle ~ 60 degs at the base of the radio jet in M87, in line with the wide outflow geometries 
found in our solutions. On the other hand, other geometries have been proposed for other 
classes of systems (e.g. Elvis 2000). 

The motivation to understand the nature and properties of outflows from accreting black 
holes goes beyond what we have just outlined, as outflows may be an efficient process by 
which black holes influence their environments. Recently, this general point has been raised 
in the speciflc context of radiatively-inefficient accretion flows by several groups. Soria et 
al. (2006a,b; see also Di Matteo, Carilli & Fabian 2001) invoked slow outflows to effect a 
feedback on the rate of Bondi accretion in several low-luminosity AGN. The wide outflow 
geometries in our solutions may suggest that the conceptual use of a simple Bondi spherical 
inflow geometry to estimate gas capture rates is inappropriate if indeed a large fraction of 
the black hole capture sphere is the site of gas outflow. The work of Nagar et al. (2005) and 
Allen et al. (2006) further suggests that feedback due to outflows associated with radiatively- 
inefficient accretion flows can reach far beyond the black hole immediate vicinity and affect 
star formation in the host galaxy, or even the surrounding hot gas when the host is located 
at the center of a group or a cluster of galaxies. It will be interesting to develop these 
arguments, and more detailed feedback scenarios (e.g. Vernaleo & Reynolds 2006), in light 
of what we have learned in our study of thermal outflows. 



6. Discussion and Conclusion 

The possibility that tangled magnetic flelds strongly limit the efficiency of thermal 
conduction in a hot flow constitutes a major theoretical uncertainty for our models. In 
recent years, a similar issue has been discussed in the context of galaxy cluster cooling 
flows. Several studies have concluded that magnetic flelds may not limit too strongly the 
efficiency of thermal conduction (e.g. Narayan & Medvedev 2001; Gruzinov 2002; see also 
Chandran & Cowley 1998; Chandran et al. 1999), but this remains somewhat of an open 
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issue. If tangled magnetic fields were to limit the effective mean free paths in the hot flow 
to Ics ^ R, conduction would likely still have important consequences for the flow structure 
and properties. If, on the other hand, tangled magnetic fields act to make l^s <^ -R, a much 
weaker role for conduction (if any) is to be expected. Explicit numerical simulations with 
anisotropic conduction along field lines in an MHD turbulent medium could greatly help 
in settling this issue (see, e.g., Cho & Lazarian 2004; Sharma et al. 2006). It should be 
emphasized that the energy transport mechanism which permits the existence of outflows 
in our solutions does not have to be thermal conduction. Any form of heat transport (e.g. 
turbulent, as emphasized by Balbus 2004) acting like thermal conduction would have similar 
consequences for the flow. 

One of several limitations of our self-similar hot accretion flow solutions with conduc- 
tion is their 1-temperature structure. Narayan & Yi (1995b) have shown how crucial it is 
to account for the 2-temperature structure of the hot flow to obtain solutions with realistic 
cooling properties. In such solutions, the expected decoupling of electron and ion temper- 
atures in the inner regions of the hot flow would modify the role of conduction. For one 
thing, ions may be partially insulated from electron conduction. In addition, self-similarity 
is broken in the 2-temperature regime and the electron temperature profile flattens to a 
sub-virial slope (e.g. Narayan et al. 1998b). In global solutions, conduction is expected to 
flatten both the polar and the radial temperature profiles (while it only affects the polar 
profile in the radially self-similar solutions discussed here). Furthermore, if a saturated form 
of conduction is appropriate, the saturated flux should be reduced in regions of the flow 
where the electrons become relativistic, since their "free streaming" velocity then becomes 
limited by the speed of light. Finally, the break-down of self-similarity implies that the radial 
contribution to conduction, which acts only to heat up the flow in our solutions, will have 
to be done at the expense of inner regions that will now be cooling by conduction. Clearly, 
it will be interesting to see how all these effects combine with each other in global models of 
hot accretion flows with conduction. 

A potentially more fundamental limitation of our work may be the use of simple fluid 
equations to describe a weakly-collisional plasma obeying a more complex set of equations. 
For instance, Balbus (2001) has argued that the dynamical structure of hot flows could be 
affected by the anisotropic character of conduction in the presence of magnetic flelds (see also 
Parrish & Stone 2005), while Quataert et al. (2002) have emphasized additional effects due 
to the anisotropic pressure tensor present in the weakly-collisional regime (see also Sharma 
et al. 2006). On the other hand, it is possible that a scale analysis could justify the use of 
simpler fluid equations to describe the dynamics of the weakly-collisional magnetized plasma 
on the largest scales of interest in the context of hot accretion flows. Such an analysis would 
clearly be valuable since essentially all the literature on hot accretion flows relies on collisional 
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fluid equations (magnetized or not). 

In recent years, much emphasis has been put on reducing the rate at which gas is accreted 
by a central black hole, within a hot flow, either by arguing that a negative Bernoulli constant 
is required for an accretion solution with imposed winds to be physically viable (Blandford 
& Begelman 1999) or by invoking a gradual reduction in the rate of mass accretion in the 
flow due to the effects of convection (Narayan et al. 2000; Quataert & Gruzinov 2000a). Our 
results, combined with the possibility that conduction may be important in hot accretion 
flows, could pose a challenge to both of these scenarios. We have already seen that the 
sign of the Bernoulli constant does not determine the existence of thermal outflows, or even 
discriminate between inflowing and outflowing regions, in our solutions with conduction. 
Our solutions also show that a radially varying m is not necessarily associated with the 
existence of such outflows. In addition, if conduction is important on all scales of interest 
in the hot flow, convection may be suppressed altogether when a displaced fluid element 
efficiently reaches thermal equilibrium with its environment through conduction. 

Obviously, our results do not exclude the existence of other classes of spontaneously 
outflowing solutions with radially varying m, density proflles deviating from the canonical 
^-3/2^ more realistic flow patterns allowing for non-zero latitudinal velocities, polar funnels, 
or a role for other radial and latitudinal energy/momentum transport mechanisms (see, e.g., 
the issues of adiabatic stability discussed in Blandford & Begelman 2004, but also Balbus 
2001 for the diabatic case with anisotropic conduction and magnetic flelds). It would be 
particularly interesting to contrast the study of these generalized solutions with the results 
for averaged flow properties emerging from simple hydro dynamical and MHD numerical 
simulations (e.g. McKinney & Gammie 2002; Proga 2005), in an effort to make the best use 
of these two complementarity approaches to the black hole accretion problem. 

If conduction is generally important in hot accretion flows around black holes, it could 
have additional interesting consequences besides the ones illustrated in our speciflc solutions. 
For instance, if the hot flow residing inside the capture radius around a black hole is able 
to heat up the ambient gas outside of that capture radius, the Bondi rate of gas capture 
(oc c~^) may be effectively reduced. While Gruzinov (1998) invoked "turbulent" heat flows 
to achieve this, the weakly-collisional conditions discussed in §2 for several galactic nuclei 
suggest that this could be done by "microscopic" conduction itself. More detailed models are 
therefore needed to address various interesting issues related to conduction, to understand 
better how conduction may affect the structure and properties of hot accretion flows and to 
determine how these effects may relate to the observational characteristics of dim accreting 
black holes. 



-21- 



We are grateful to S. Balbus, A. Beloborodov, W. Borland, J. -P. Lasota, R. Narayan, 
E. Quataert and E. Spiegel for useful discussions and comments. We also thank M. A. Morse 
for assistance in producing the mass flux color flgure. This research was supported in part 
by the National Science Foundation under Grant No. PHY99-07949 (at KITP). 

APPENDICES 



A. Prescriptions for Conduction 

In §3, we have seen that a saturated form of conduction scales in such a way that 
it permits ID radially self-similar solutions. Consequently, we flrst attempted to find 2D 
radially self-similar solutions with a saturated conduction flux written as 
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which depends only on the direction of the local temperature gradient, not its magnitude. 
From taking minus the divergence, we obtained the following conduction term in the energy 
equation: 
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within the framework of our self-similar solutions. 



We have found this prescription to be numerically unstable, or at least difficult to work 
with. This may be due to the explicit coupling between the radial and polar conduction 
terms that appears within this prescription (which does not seem well suited to our separable, 
radially self-similar solutions). It may also be caused by the structure of the mathematical 
operator resulting from this prescription: it no longer acts as a diffusion process. 

Consequently, in our 2D treatment (in §4), we returned to a standard diffusive operator 
for conduction, with a conductivity coefficient scaled with radius to preserve self-similarity. 
This results in the following conduction term in the energy equation (see also Eq. [23]): 
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From a comparison of these two prescriptions, it is clear that, in the hmit of strong 
conduction, when latitudinal profiles of most flow quantities become nearly fiat, the two 
conduction parameters, $s and Aq, are simply related to each by Aq ~ bp{6)cs{6) $s- Using 
the simple scalings of Appendix B, the relation becomes 



Ao^ ^^$s = 8.4$„ (A4) 

where we have used m = 1, a = 0.1 for the last equality. This correspondence is helpful in 
showing that, in all our 2D solutions, the non-rotating limit is reached for $s equivalent values 
which are <^ 1, in good agreement with our ID results (in §3). In other words, even in 2D 
solutions, a level of conduction well below that of the (non-relativistic) saturated conduction 
flux is sufficient to alter the structure and properties of hot accretion flows substantially. 



B. Maximum Value of the Conduction Coefficient 

As we have seen in §4, for a given value of the flow parameter e', there is a maximum value 
of the conduction parameter, Aq, at which the solutions reach the non- rotating, Bondi-like 
limit. No solutions of physical significance exist beyond this maximum value since they would 
require imaginary values of Vt, as we shall now establish. Let us evaluate Eqs. (20-23), at 
the fiow equator, using adequate boundary conditions, in the limit of low viscosity {av <^ 1), 
small radial velocity (f ^ ^ 1) and under the assumption that second order derivatives with 
respect to 9 become small (which is generally accurate for solutions approaching the non- 
rotating limit). We obtain: 

(Bl) 
(B2) 

(B3) 

The first equation establishes that, in the strong conduction limit, Cg ~ |. The second 
equation yields, in combination with the first, Vr ~ ""f'^? ^^ the same limit. To be exact, 
this is not the limit f] ^ 0, but rather (PQ/dO"^ ~ 0. The actual non-rotating limit occurs at 
a slightly lower value of Vr ^ "fct- The third equation places the explicit constraint on Aq 
for the non-rotating limit, since the left-hand side must be positive for the value of Q to be 
real. This fixes the maximum value of the conduction parameter Aq to 
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where we have used the crude approximation that momentum (pvr) profiles are uniform 
along 6, which results in p{6) ~ 5m/(67rQ;). This estimate of the maximum value of Aq 
agrees well with our numerical results (e.g. in Fig. 6). 



C. Comparison with Anderson (1987) 

Our results on 2D hot accretion flows with conduction share strong similarities with 
those of Anderson (1987). First, let us mention that we have succeeded in reproducing many 
of the results reported by this author, using the same prescription for conduction ( i.e. with a 
conduction coefficient which scales not only with radius, but also with the local flow viscosity 
and density). With this prescription, the conduction term in the energy equation becomes: 
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Qa = aAi 



^(*'<<*' + 5^5? =-»^(»> ,9 



(ci; 



Returning to our analysis, there are still differences between our work and that of 
Anderson (1987), both in the formalism and in the results. Anderson works in the limit of 
"weak shear viscosity," with several viscous terms dropped from the dynamical equations 
(as compared to ours). He adopts a polytropic equation of state, a prescription for viscosity 
which differs from Shakura-Sunyaev and there is no inertial term in his radial momentum 
equation. As a result, his set of differential equations is reduced to a third order system, as 
opposed to a seventh order one in our case. Despite these technical differences, our results 
are qualitatively consistent with many of his conclusions. Differences appear mostly at the 
quantitative level, for instance in the rates of mass outflow (which appear a few times larger 
in our solutions). One qualitative feature which has apparently been missed by Anderson 
(1987) is the existence of solutions with conical outflow geometries in sub- regions of the 
parameter space. 

From this comparison, we can safely conclude that many of the qualitative features 
exhibited by our solutions should be relatively robust to changes in the various prescriptions 
adopted (e.g. conduction and viscosity). 
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Table 1: Observational Constraints on Collisionality in Nearby Galactic Nuclei 
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Fig. 7. — Latitudinal profiles of the advection (a), conduction (b) and viscous (c) terms in 
the 2D self-similar energy equation, for five solutions with a = 0.1, e' = 10 and varying 
degrees of conduction (as measured by Ag). Note the different vertical scales. 
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Fig. 10. — Latitudinal profiles of the dimensionless Bernoulli parameter, h{9) (solid lines), 
and of the radial velocity, Vr{6)/vK (multiplied by an arbitrary factor 5; dashed lines) in 2D 
self-similar solutions with 7 = 1.5, / = 0.033 (hence e' ~ 10) and two values of the scaled 
conduction parameter, Ap/e'. The two panels show solutions of special interest, with a conical 
outflow geometry (right panel) and with a marginal outflowing region having b{6) < (left 
panel). Horizontal dotted lines were added in each panel to help identify inflowing vs. 
outflowing regions. 



